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I.  Introduction 

The  recent  discovery  that  submonolayer  films  of  adsorbates  on  some 
substrates  have  two-dimensional  phases  has  opened  a  new  field:  the  study  of 
two-dimensional  matter,  its  phases,  and  the  transitions  between  them.  A 
common  way  to  realize  two-dimensional  systems  is  to  physisorb  gaseous 
particles  (atoms  or  simple  molecules)  onto  solid  surfaces  (substrates). 

Dependent  on  particle  and  substrates  properties,  the  adsorbate  condenses  as 
a  homogeneous  film  which  might  reveal  two-dimensional  properties.  Often  the 
only  task  the  substrate  is  expected  to  accomplish  is  to  force  the  adsorbed 
particles  into  a  plane,  so  that  the  substrate  disturbs  the  adsorbed  system 
as  little  as  possible.  In  order  to  investigate  experimentally  the  physical 
properties  of  two-dimensional  systems,  one  must  have  substrates  of  high 
specific  surface  areas  such  as  graphite  (and  its  modifications).  The 
essential  task  in  a  thermodynamical  investigation  of  a  physisorbed  two- 
dimensional  system  is  the  construction  of  a  complete  group  of  adsorption 
isotherms.  In  practice,  for  an  accurate  determination  of  the  entropy, 
calorimetric  experiments  to  determine  the  specific-heat  signatures  are 
performed.  Such  experiments  offer  the  possibility  of  studying  low¬ 
dimensional  order-disorder  phenomena  such  as  structural  and  melting  phase 
transitions. 

With  respect  to  the  order  of  a  melting  transition,  the  situation 
depends  strongly  on  whether  one  considers  three-  or  two-dimensional  melting. 

In  the  former  case,  there  is  no  doubt  that  melting  is  always  a  discontinuous 
event  at  which  the  system  absorbs  latent  heat  from  the  surroundings.  In  the 
two-dimensional  case,  the  question  as  to  whether  melting  is  a  discontinuous 
or  continuous  process  is  still  open. 


In  the  present  paper  we  use  the  McQuistan-Hock  model  to  compute  heat- 

capacity  signatures  for  two-dimensional  systems  of  rare  gases  adsorbed  on 

graphite.  Experimental  studies  on  these  systems  have  demonstrated  the 
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existence  of  ordering  and  phase  transitions  and  have  identified  regions  of 

two-dimensional  liquid-vapor  as  well  as  liquid-solid  coexistence, 

3 

accompanied  by  a  two-dimensional  critical  point.  In  the  next  section, 
after  a  brief  overview  of  the  McQuistan-Hock  model,  we  present  our 
calculations  and  results.  The  third  and  last  section  is  devoted  to  our 
conclusions. 

II.  Theory  and  Calculations 

In  a  recent  paper, *  McQuistan  and  Hock  developed  an  exact  solution  for 
the  distribution  function  of  q  indistinguishable  particles  on  a  2  *  N 
lattice,  where  N  is  a  positive  integer  number.  The  grand-canonical 
partition  function  is  written  as 


Zjj(x,y,z,t)  »  £  FqN(x,y,t)s 
q-0 


where 
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in  which  p  is  the  chemical  potential  of  the  adsorbed  particles,  p  =  (k_T)  * 
with  T  the  absolute  temperature  and  kg  the  Boltzmann  constant,  Vqq  is  the 
interaction  between  two  nearest-neighbor  vacant  sites,  is  the 
interaction  between  two  nearest-neighbor  adparticles,  and  Vq  is  the 
interaction  between  a  particle  and  the  surface.  The  number  of  unique  ways  q 
indistinguishable  particles  can  be  arranged  on  a  rectangular  2  *  N  lattice 
to  form  n^  occupied  nearest-neighbor  pairs  and  n^Q  vacant  nearest-neighbor 
pairs  is  given  by  A[N,q,ngQ,n^]  in  Eq.  (2)  by  making  use  of  a  15-term 
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recursion  relation.  In  a  subsequent  paper,  Hock  and  McQuistan  computed  the 
adsorption  isotherms  for  what  we  might  call  the  "rigid"  lattice  in  which 
there  is  no  interaction  between  two  nearest-neighbor  vacant  sites:  s  0. 

They  showed  that  the  coverage  as  a  function  of  the  spreading  pressure  does 
not  exhibit  a  first-order  phase  transition. 

The  McQuistan-Hock  model  is  very  interesting  because  of  its  simplicity, 
because  it  takes  into  account  various  interaction  terms  (Vq>voq  and  V^), 
and  because  it  is  exactly  solvable.  We  will  show  here  that  the  model  is 
applicable  to  realistic  systems  by  comparing  the  critical  temperatures 
obtained  from  it  with  the  critical  temperatures  measured  in  systems  such  as 
rare  gases  adsorbed  on  graphite. 

In  the  McQuistan-Hock  model,  the  grand  canonical  partition  function 
[Eq.  (1)],  as  N  ■»  «,  is  given  by 

Z(N,p,p)  *  »  (7) 

where 


rCnj) 

ki  '  -  F7^>  • 


(8) 


and  is  the  smallest  root  of 


*'l 


iffl 
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'•t; 


s(n)  =  £  bj  ■  0  . 

j-0 


In  Eq.  (8), 


:<n)  -[  VJ 
j-0 


s'(n)  = 


(— ) 

on  x»y>z 


The  coefficients  aj(x,y,z)  and  bj(x,y,z)  in  Eqs.  (9)  and  (10)  are  given  in 
Eqs.  (9)  and  (10)  of  Ref.  4.  The  spreading  pressure  p  can  be  obtained  as 


1  .  20s' 

*  ■  H  1“"ir 


where  6  is  the  coverage  as  N  •*  •»,  and 


s(h)  "  (ff),  v  „  * 

3z  x,y,n 


In  order  to  study  the  behavior  of  the  heat  capacity  as  a  function  of 


the  temperature,  keeping  the  spreading  pressure  p  and  the  number  of 


adparticles  q  as  constant  it  is  most  convenient  to  work  within  the 


isothermal- isobaric  ensemble,  whose  partition  function  will  be  denoted  by  A 


*  A(q,0,p).  The  Gibbs  free  energy,  entropy  and  heat  capacity  are  given  by 


G(q,P,p)  - 


£ln*  • 


S(q,P,p) 


.  Ml 


C(q,P,p) 


V2  (f  '  (f)2>  ’ 


where  C(q,g,p)  has  a  maximum  at  the  critical  temperature  (0  *  g  ),  i.e., 
the  solution  of  the  equation 

2A(A')2  +  3gAA*  A"  -  2g(A’)3  -  gA2AH  -  2A2Am  -  0  .  (17) 


In  Eqs.  (15)-(17),  the  prime  denotes  the  derivative  with  respect  to  g.  In 
order  to  find  A(q,g,p),  we  start  from  its  definition 

A(q,g,p)  s  Aq(g,p)  3  ^  Fq(H,g)nN 

N-l 

00 

-  tq  £  fq(N,g)n N  ,  (18) 

N-l 

where  n  *  e  2^*\  Multiplying  by  z^,  summing  over  q  from  0  to  2N  and  using 
Eqs.  (3)  and  (4)  given  in  Ref.  4,  the  following  recursion  relation  is 
obtained  for  Aq(g,p): 


Vq  *  ’  V,-l  '  Vq-2  '  d3‘q-3  +  c36q,3  • 

(19) 

where  6  -  is  the  Kronecker  delta  function, 

<!*■» 

d0  -  l  -  y3n  . 

(20a) 

dj  -  -n(xydQ  +  y3n  +  l)t  , 

(20b) 

d2  -  -n(*3d0  +  xyn)t2  , 

(20c) 

d3  -  xn2(xy-l)lx2dQ  +  yn(2xy-l)]t3  , 

(20d) 

and  c^  is  given  by  Eq.  (24d)  below.  From  Eq.  (18)  and  from  the  recursion 
relation  for  fq(N,g)  as  given  in  Ref.  4,  we  can  compute  the  initial 
condition  for  Aq(g,p).  We  obtain,  provided  Vqq  >  -2p/3,  that 
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a1  ■  2n  +  4ynAQ  +  2y2n^Q  , 

A2  *  xn(l+2yn)  +  2nAQ(2xy2n  +  y2n  +  X  +  1/y) 

+  nA^(2xy3n  +  4y3n  +  x  +  4)  +  2ynA3(l  +  y3r)) 


(21b) 


(21c) 


Let  us  define 


\  .  u(-z) 

q  v(z) 


where 


u(z)  -  £  Cj  z^ 
j-0 

3 

v(z)  »  ^  dj  z^  . 
j=0 


(23a) 


(23b) 


Here  the  dj's  are  given  by  Eqs.  (20)  and  the  Cj's  are  given  by 


cn  -  yn 


Cj  -  nt  [2dQ  +  yn(4y  -  xy  -1)] 
c2  ■  nt2[xd0  +  xyn(2-x2)] 

c3  -  xn2t3{[x2yn  +  4yq(l-x)](xy-l)  +  (4-y-2x)(x-xy3n+y2q) 


If  we  rewrite  Eq.  (22)  as 


(24a) 

(24b) 

(24c) 

-  1)  .  (24d) 


u(z)  »  v  . 

v(z)  "  K  + 


Cq  +  CjZ  +  C2z 

1  +  Dxz  +  D2z2  +  D3z3 


■  K  +  z 


with 


K  *  *Vd0  » 


d0  Cj  "  cj+l  d0  '  C0  dj+l 


we  obtain 


3 

Aq  "  I  x(zj)  zjq  *  0  <  q  -  2N 


where  the  zj's  are  the  roots  of  the  equation 


6(z) 


(29) 


(30) 


and 


X(z)  =  -Y(z)/6*(z)  .  (31) 

Here  the  prime  denotes  the  derivative  with  respect  to  z. 

It  is  worth  noting  that  for  the  special  case  in  which  x  *  y  »  1,  = 

-  0  and  the  solution  of  Eq.  (29)  gives 


z 


± 


~n_*.Vn 

nt 


(32) 


and  the  postive  root  is  always  smaller  than  the  absolute  value  of  the 
negative  root,  so  that  for  large  q, 


A 

q 


(Wo)q 

2(1-Vn)q+1  ’ 


and  the  heat  capacity  per  adparticle  is 


kBa 

sinh^a 


(33) 


(34) 


y 

where  a  =  flp/2.  Vie  note  that  this  is  the  same  result  that  would  have  been 
obtained  for  a  system  of  independent,  distinguishable  particles,  each  of 
which  has  only  two  accessible  states:  the  particle  is  on  the  site  or  it  is 
not. 

For  the  more  general  case  in  which  *  0,  we  have  solved  Eq.  (30) 
numerically  choosing  values  for  the  interaction  constants  that  are  likely  to 
model  rare-gas/graphite  systems.  In  Table  I,  the  chosen  values  are 
displayed  together  with  the  literature  reference  from  where  they  have  been 
taken.  Several  studies  of  the  Ne,  Ar,  Xe/graphite  systems  have  reported  a 
broad  anomaly  in  the  specific  heat  at  temperature  values  given  in  Table  I, 
namely  at  15.7  K  for  Ne,  55  K  for  Ar  and  118  K  for  Xe.  In  Fig.  1  we  present 
the  specific  heat  as  a  function  of  temperature  for  values  of  the  coverage  of 
0.1,  0.3  and  0.5.  We  see  that  maxima  in  the  specific  heat  are  observed  in 
the  range  of  the  experimental  critical  temperatures. 

III.  Discussion  and  Conclusions 

From  Fig.  1  we  can  see  that  maxima  in  the  heat-capacity  signatures  are 
found  around  10-15  K  for  Ne,  40-52  K  for  Ar  and  80-102  K  for  Xe.  Let  us 
compare  these  results  with  some  experimental  findings  and  with  some  other 
model  c-^culations. 

g 

Steele  and  Karl  reported  peaks  in  the  specific-heat  measurements  of  Ne 
adsorbed  on  graphitized  carbon  black  powder,  where  at  half  coverage  they 
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found  a  peak  at  16.1  K.  Antoniou  in  turn  found  a  peak  at  12  K.  Huff  and 
Dash***  made  the  first  set  of  measurement  of  Ne  adsorbed  on  Grafoil,  for 
submonolayer  coverages  in  the  temperature  range  of  2-20  K  and  found  peaks 
between  12  and  15  K.  The  anomalies  they  obtained  have  been  confirmed  by 
Rapp  et  al,**  who  found,  besides  a  first  peak  at  13.6  K,  a  second  one  at 


around  16  K  which  moves  up  in  temperature  as  the  coverage  increases,  and 
then  disappears  at  higher  coverages. 

12  13 

Recent  experimental  studies  of  neutron  scattering,  *  heat 

14  15 

capacities  and  synchrotron  x-ray  scattering  on  submonolayers  of  Ar  on 

graphite  were  all  interpreted  as  being  consistent  with  a  continuous  melting 

process.  Migone  et  al^  reported  heat-capacity  measurements  of  submonolayer 

Ar  on  Graphite  foam,  with  anomalies  at  47.2,  49.5  and  55  K,  and  interpreted 

the  first  peak  as  a  "weak"  first-order  transition.  For  Xe  the  situation  is 

similar:  melting  of  Xe  physisorbed  on  graphite  has  been  interpreted  both  as 

17  18  19 

continuous  ’  and  first-order  melting. 

The  interesting  aspect  of  the  McQuistan-Hock  model  is  that  it  does  not 
allow  for  any  first-order  transition.  Yet,  when  applied  to  "realistic" 
systems,  it  predicts  maxima  in  the  heat-capacity  signatures  at  temperatures 
very  close  to  the  experimental  ones.  A  single  model  that  forbids  a  first- 
order  transition  but  displays  maxima  in  the  heat-capacity  curve  is  the  one¬ 
dimensional  lattice  with  nearest-neighbor  interactions.  Unfortunately,  this 
model  predicts  those  maxima  at  temperature  values  well  below  the 
experimental  data.  A  two-dimensional  model  such  as  the  square  lattice 
predicts  second-order  critical  temperatures  above  the  experimental  data, 
while  the  hexagonal  lattice  gives  better  agreement  with  experiments. 

In  Table  II  we  report  the  ratio  R  =  as  Predicted  by  various 

models,  as  well  as  the  experimental  average  value  K  of  this  ratio  as 

R.tr 

obtained  from  Table  I  ( l"p~l  *  10%).  The  predictions  from  the  McQuistan- 
Hock  model  are  also  reported  as  averages  because,  as  can  be  seen  from  Fig. 

1,  R  is  not  exactly  constant.  From  this  table  it  is  interesting  to  see  that 
the  McQuistan-Hock  model  predicts  values  of  Tc  that  are  in  good  agreement 


with  the  experimental  measurements.  The  model  does  not  allow  for  first- 
order  transitions,  nor  does  the  linear-lattice  model  with  nearest-neighbor 
interactions.  On  the  other  hand,  the  linear-lattice  model  predicts  maxima 
in  disagreement  with  experiments.  In  good  agreement  with  experiments  are 
also  the  anomalies  predicted  by  the  square- lattice  and  by  the  hexagonal 
models,  the  difference  between  the  two  being  what  we  call  the  coordination 
number  (CN  in  Table  II),  i.e.,  the  number  of  nearest  neighbors  attached  to 
each  site:  '  the  experimental  Tc's  lie  in  between  the  predictions  from  the 
models  with  CN  ■  3  and  CN  •  4. 

Our  conclusions  are  the  following:  (i)  The  McQuistan-Hock  model  has 
the  advantage  of  being  exactly  solvable  and  of  giving  good  agreement  between 
computed  and  experimental  temperatures  at  which  the  maxima  in  the  heat- 
capacity  signatures  occur,  (ii)  The  model  contains  very  few  parameters, 
usually  Vq,  Vqq  and  V^j  yet  those  parameters  are  sufficient  to  determine 
the  "gross"  features  of  the  statistics  of  adsorption.  This  implies  that  the 
details  of  the  various  interactions  occurring  in  real  systems  do  not  seem  to 
play  an  important  role,  (iii)  The  model  contains  a  parameter,  Vqq,  which 
accounts  for  distortions  of  the  substrate,  an  effect  that  has  been  invoked 
by  various  authors^’ ^  to  explain  some  details  of  the  experimental  heat- 
capacity  signatures.  It  would  seem  worthwhile  to  explore  the  dependence  of 
the  adsorption  thermodynamics  on  Vqq.  (iv)  The  impossibility  of  showing  a 
first-order  transition,  on  one  hand,  and  the  good  agreement  between  the 
computed  and  measured  critical  temperatures  suggests  that  the  order  of  the 
transition  seen  at  temperatures  around  the  Tc-values  of  Table  I  might  not 
necessarily  be  of  first  order,  (v)  An  important  role  seems  to  be  played  by 
the  number  of  nearest  neighbors  attached  to  each  site:  the  McQuistan-Hock 
model  (with  CN  ■  3)  represents  a  logical  extension  from  a  one-dimensional  to 


the  complete  two-dimensional  system.  It  would  therefore  be  interesting  to 
examine  the  influence  of  the  average  number  of  nearest  neighbors  on  the  heat 
capacity  as  a  function  of  temperature. 
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Table  I.  Adparticle-surface  (Vg)  and  adparticle-adparticle  (V^) 
interactions  in  Kelvin.  Tc  is  the  critical  temperature. 


*vo 

b  V 

11 

c  T 

c 

d  T 

c 

Ne 

-378 

-34.6 

15.7 

14.5 

Ar 

-1103 

-120 

55 

51.0 

Xe 

-1912 

-236 

118 

102.0 

a  From  Ref. 
^  From  Ref. 

Q 

From  Ref. 
d  From  this 


5 

6 
7 

study  at  coverage  value  of  0.5 


Table  II.  Values  of  the  ratio  R  *  |Vjj/Tc|  from  various  models  and  from 
experiments  (EXP);  ID  •  linear  lattice,  HL  ■  hexagonal  lattice, 
MQH  *  McQuistan-Hock  lattice,  SL  ■  square  lattice.  For  EXP  and 
MQH  the  average  values  are  reported.  The  coordination  number  is 
indicated  in  parentheses. 

Systems  R 


ID  (CN  -  2) 

4.80 

HL  (CN  -  3) 

2.63 

MQH  (CN  »  3) 

2.35 

EXP 

2.13 

SL  (CN  -4) 

1.76 

.*  ■»  Vrv:v,v  y  v  v„ 
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Figure  Caption 

1.  Specific  heat  for  Ne,  Ar  and  Xe  adsorbed  on  graphite  using  the 

interaction  parameters  given  in  Table  I  and  the  coverage  values  of  0.1, 
0.3  and  0.5. 
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